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1.  INTRODUCTION 


It  is  a  common  practice  in  metallurgical  fatigue  testing  that  m  (  >1) 
levels  of  stress  are  assumed  to  be  employed  with  n  specimens  tested  at 
each  of  m  levels  of  stress.  The  underlying  distribution  of  stress  cycles 
to  failure  at  each  level  (i*"*1  level)  is  assumed  to  be  a  two-parameter 
Weibull  distribution. 


Fi(x)  =  1  -  exp  {-(x/bi)Cil  ,  x  >  0,  (1) 

where  c^  and  b^  are  the  shape  and  scale  parameters  respectively. 

From  the  fatigue  data  x.,,  x.„ .  x.  (under  ith  level  of 

ll  1 2  in 

stress) ,  the  unknown  parameters  c^  and  b^  are  usually  estimated  by  the 

maximum  likelihood  method.  The  maximum  likelihood  estimators  c  and  b 

i  l 

are  found  to  be  solutions  of  the  following  equations; 

n 


n 

7  x . .  i  An  x . .  .  7  An  x . . 

j=i  iJ  1  _  j£i _ ^ 


n 


=  0 


l  x..ci 

U  IT 


j  =  l  1J 


c  . 
1 


n 


and 


/s  n  1/c. 

b.  =  {l/n7x..Ci) 

1  j=i 


(2) 


(3) 


Cohen  [1]  has  suggested  a  method  of  numerical  computation  of  c.  and 
Thoman,  Bain  and  Antle  [5]  have  given  statistical  inferences  on  c  and  b  . 

l  l 

However  there  is  a  conjecture  and  strong  experimental  evidences  that 

the  shape  parameters  are  independent  of  applied  level  of  stress  (see 

Hahn  and  Kim  [2]  and  Lipson,  She.th  and  Desney  [3]).  This  implies  then 

that  c,  =  c„  =  . . .  =  c  =  c  and  it  is  needed  to  estimate  the  common 
12  m 

shape  parameter  c  on  the  basis  of  the  pooled  data  x  ,  i  =  l,2,...,m  and 
j  =  l,2,...,n.  Statistical  tests  for  the  equality  of  the  shape  parameters 
c^  were  given  in  Thoman  and  Bain  [4]  and  Schafer  and  Sheffied  [6]  . 

Wolff  and  Lemon  [7]  have  recently  considered  various  pooled  estimation 
methods  of  c  to  analyze  composite  materials  data,  and  for  the  case  of  m  =  2, 
pooling  techniques  were  involed  in  the  study  of  two-sample  tests  by  Schafer 
and  Sheffied  [6]  and  Thoman  and  Bain  [4]. 


1 


The  problems  of  estimations  and  testing  hypotheses  regarding  to  the 
common  shape  parameter  c  and  scale  parameters  in  the  Weibull  distribu¬ 
tions  are  considered  in  this  report.  The  following  results  are  obtained: 

1.  Three  methods  of  pooled  estimations  of  c  (Averaging  M.L.  estimation, 
M.L.E.  by  normalization  and  Joint  M.L.  estimation)  are  given. 

These  methods  were  introduced  by  Wolff  and  Lemon  [7]  but  included 
here  for  completeness. 

2.  Exact  confidence  intervals  and  testing  hypotheses  for  c,  depending  on 
the  methods  of  pooled  estimations  of  c,  are  presented. 

3.  Comparison  of  the  bias  factors  of  the  three  pooled  estimators  of  c 
is  given. 

4.  Exact  confidence  intervals  and  testing  hypotheses  for  b^  (scale 
parameter  at  each  level  of  stress) ,  depending  on  the  method  of  pooled 
estimators  of  c,  are  presented. 

2.  POOLED  ESTIMATIONS  OF  c 


Let  x.,,  x._,  ...,  x.  be  a  random  sample  of  size  n  taken  from  a 
il  i2  in 

Weibull  distribution  with  the  shape  parameter  c  and  scale  parameter  b^, 
for  i  =  l,2,...,m.  Three  methods  of  pooled  estimation  of  c  are  presented 
in  this  section. 


(A) 


Averaging  M.L.E.  c  :  Let  c^  be  the  M.L.E.  of  c^  based  on  the 

observations  x.,,  x.„,  ...,  x.  .  Since  the  averaging  procedure  is 
li  lZ  in 

a  common  technique  in  estimations,  the  averaging  M.L.  estimator  c 
of  c  is  defined  by 


m 

c  =  1/m  I  c.  .  (4) 

i=l  1 


(B)  M.L.E.  by  normalization  c  :  Since  the  shape  parameter  c  is  free 

from  scale  changes  (normalization),  pooled  estimation  of  c  can  be 

obtained  by  normalizing  the  data,i.e.  by  letting  y  =  x../b.,  for 

i  J  i  3  i'v, 

i  =  l,2,...,m  and  j  =  l,2,...,n.  The  M.L.E.  by  normalization  c  is  a 
solution  of  the  following  M.L.  equation  for  the  normalized  data  : 


2 


m  n 


it 


£  y  .  ,L  5-n  y  .  ,  ,  £  £  £n  y .  . 

i=l  1=1  1J  1J  1  i=l  1=1  1J 


m  n 


l  l 


% 

c 


=  0 


(5) 


m  n 


i=l  j=l 


ij 


(C) 


When  normalizing  the  data,  the  scale  parameters  b.  are  usually  not 

^  A 

known,  and  hence  they  may  be  replaced  by  the  M.L.E.  b^  based  on 

the  observations  x._,  x, x.  . 

ll  i2  in 

Joint  M.L.  estimator  c  :  The  M.L.  equation  for  pooled  data 
x_  ,  i  =  1,2,. . .  ,m  and  j  =  1,2,. .  .  ,n  ,  is 


V  c  n 

,  /  X.  .  to  X,. 

>1  *J  lJ 


iii  L 


m  n 

r  r 

L 


X  . 


—  _  1-1  j-1 

A 

c  n 


/  fcn  x.  . 
L  ij 


=  0  , 


(6) 


j  =  l 


ij 


and  a  solution  c  of  (6)  is  called  the  joint  M.L.  estimator  of  c. 


For  notational  convenience,  let  c  be  representing  any  one  of  pooled 
estimators  of  c  given  in  (A),  (B)  and  (C).  We  also  denote  the  corres¬ 


ponding  M.L.E.  of  b^  by 


b .  =  [  1/n  y  x . . 


*  ,  * 
c  .  1/c 


(7) 


j-1 


for  i  =  1,2, ... ,m  . 


3.  CONFIDENCE  INTERVALS  FOR  c 

In  what  follows,  and  are  used  to  denote  the  pooled  esti¬ 

mation  of  c,  given  in  (A) ,  (B)  and  (C) ,  when  in  fact  the  sample  is  from 
a  Weibull  distribution  with  b  =  1  and  c  =  1,  i.e.  a  standard  exponential 

•Jf  _  ^  A 

distribution.  Let  c^  be  representing  any  one  of  c^  and  c 

The  key  result  for  the  estimation  of  c  is  the  following  theorem, 
which  can  be  proved  by  using  the  same  method  given  in  the  proof  of 
Theorem  A  of  Thoman,  Bain  and  Antle  [5].  Hence  the  proof  of  the  following 
theorem  is  omitted  . 
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THEOREM  1.  c*/c  is  distributed  independently  of  b.  ,  b„ . b 

12  m 

and  c  and  has  the  same  distribution  as 

The  distributions  of  c^,  and  c  ^  were  obtained  by  the  Monte 

Carlo  method-  These  distributions  were  based  on  the  simulations  of 
10,000  random  samples  of  size  n  x  m  which  was  performed  at  the  Air  Force 
Aeronautical  Systems  Division  (Wri ght -Patterson  Air  Force  Base)  on  the 

—  'Xz  a 

CDC  6600.  Percentage  points  of  the  distributions  of  c^>  c  and  c^ 

are  given  in  Tables  1,  2,  and  3  respectively  for  m  =  2, 3, 4, 5, 7  and 

n  =  5,6,7,8,9,10,12,14,16,18,20,25,30,50.  These  results  can  be  used 

to  construct  confidence  intervals  for  c  when  b,  ,  b„,  . . . ,b  are 

12  m 

unknown.  100(1  -  y)  percent  confidence  intervals  will  be  of  the  form 
(  c*/£  ,  c*/£<2  )  where  £  and  £9  are  from  Tables  1,  2,  and  3  such 
that  P^_i  £^  <  c*  <  }  =  1  ~  y  . 

"k  k 

The  generated  distributions  of  c.^  provides  the  factors  B  (n) 
such  that  E[B*(n)  c*]  =  c  .  We  note  that  B*(n)  represents  B(n), 

B(n)  or  B(n),  corresponding  to  c,  c  or  c.  These  unbiasing 
factors  are  given  in  Table  4. 


The  problem  of  testing  hypotheses  Hq  :  c  =  cq  against  :  c  =  c^ 

with  the  level  of  significance  Y  ,  can  be  solved  by  using  Tables  1,  2, 

and  3.  When  c,  >  c  ,  the  distribution  of  c*/c  yields  the  critical 
1  o  o 

region  (c  £  ,  00  ) .  The  power  of  the  test  is 

o  1~Y 


P  i  c*  >  c  l, 
r  o  1-Y 


H1  }  =  Pr(  C11  "  *l-YCo/cl  }  * 


which  is  independent  of  b, ,  b„,...,  b  and  depends  only  on  c  /c, ,  y 

12m  o  1 

and  n.  A  similar  approach  can  be  taken  for  the  case  that  c,  <  c  . 

1  o 

The  powers  of  the  test  as  a  function  of  c,/c  were  obtained  and  given 

1  o 

in  Figure  1  (m  =  2,  n  =  5)  and  Figure  2  (m  =  4,  n  =  14).  It  is  noted 
that  c  gives  most  powerful  test. 


4.  CONFIDENCE  INTERVALS  FOR  b 

Let  b  represent  any  one  of  b, ,  b  ,...,  b  ,  where  b  is  the  scale 

12m  i 

parameter  of  Weibull  distribution  under  i^  level  of  stress.  This  b 

—  n,  a. 

can  be  estimated  by  b  ,  representing  any  one  of  b,  b  or  b ,  given  in  (7). 
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£  —  % 

Recall  that  c  (c,  c  or  c)  appeared  in  (7)  is  an  estimator  of  c  from 
the  pooled  data.  As  before  b*^  will  denote  the  M.L.E.  of  b,  as  given 
in  (7),  when  in  fact  the  sampling  is  from  a  Weibull  distribution  with 
b  =  1  and  c  =  1  . 

THEOREM  2 .  c*£n(b*/b)  is  independent  of  b  and  c  and  has  the  same 

distribution  as  c*^£n(b*^)  • 

The  proof  of  Theorem  2  is  omitted  since  it  is  simply  done  by 
using  the  same  approach  as  in  the  proof  of  Theorem  B  of  Thoman,  Bain 
and  Antle  [5]. 

—  —  <\j  %  ~  ~ 

The  distributions  of  >  c^JlnCb^)  and  cn^n(b  )  were 

obtained  by  Monte  Carlo  methods,  based  on  the  results  of  10,000  random 
samples  of  size  n  x  m  .  The  percentage  points  of  these  distributions 
are  given  in  Tables  5,  6,  and  7. 


100  (1-Y)  percent  confidence  interval  for  b  can  now  be  constructed 
and  will  be  of  the  form 


,  ,  *  -£0/c  ,  * 

(  b  e  2  ,  b  e 


■2  /r* 
1/C  ) 


where 


and  &2  are  obtained  from  Tables  5,  6 
Pr(  \  <  CalnO.jp  <  l2  )  ■  1 


or  7  such  that 
-  7  • 


5.  EXAMPLES 

We  consider  the  following  fatigue  life  data,  given  in  Kim  and 

Park  [8],  of  composite  material  T300/5208  graphite  .’poxy  laminate  with 

[0/90/+45]  orientation, 
s 

(1)  Under  stress  level  345MPa,  fatigue  life  in  cycle  : 


293,000 

443,870 

661,090 

923,840 

1,340,070 

364,200 

491,800 

671,540 

943,300 

1,367,890 

367,580 

539,980 

704,870 

996,170 

1,488,150 

369,890 

614,960 

764,680 

1,013,630 

1,809,060 

412,200 

631,230 

778,380 

1,104,570 

3,690,560 

429,960 

646,370 

793,340 

1,333,390 
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Under  stress  level 

414MPa,  fatigue 

life  in 

cycle  : 

7,180 

17,950 

28,440 

37,330 

52,350 

10,190 

21,270 

28,760 

37,560 

66,410 

10,300 

22,080 

31,110 

38,400 

77,130 

12,740 

22,400 

33,690 

39,480 

78,720 

15,760 

22,550 

34,970 

43,680 

101,300 

17,230 

24,570 

35,470 

47,640 

Denote  b^,  c^  and  b2 

,  c2  the  scale 

and  shape  parameters 

Weibull  fatigue  life  distributions  under  stress  levels  345MPa  and 
414MPa  respectively.  The  maximum  likelihood  estimates  of  the  above 
parameters  are  : 

b  =  1,009,350  c1  =  1.58 

b2  =  39,560  c  2  =  1.71 

In  testing  hypotheses  Hq  :  =  c2  against  ,  we  must 

accept  the  Hq  since  the  critical  region  is  c2^c\  >  =  0.1) 

according  to  Thoman  and  Bain  [4]. 

_  %  /V 

The  pooled  estimates  of  c  are  :  c  =  1.645,  c  =  1.636, and  c  =  1.638. 
(one-sided)  95  percent  confidence  intervals  for  c  are  :  (  1.315,  00 ) , 
(  1.349,  °°),and  (  1.  333,  °°)  from  Tables  1,  2,  and  3  respectively. 
The  M.L.E.  of  b^  (according  to  (7))  are  : 

b  =  1,023,575  b2  =  39,121 

fc^  =  1,021,587  $2  =  39,064 

A  A 

b1  =  1,022,028  b2  =  39,077  , 

and  (one-sided)  95  percent  confidence  intervals  for  b^  and  b2>  using 
Table  6  of  the  joint  M.L.E. ,  are  :  (  846,248,  00  )  for  b^  and 

(  32,359,  °°  )  for  b  . 
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